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Abstract 

In order to describe the dynamics of monochromatic surface waves 
in deep water, we derive a nonhnear and dispersive system of equa- 
tions for the free surface elevation and the free surface velocity from 
the Euler equations in infinite depth. Prom it, and using a multi- 
scale perturbative methods, an asymptotic model for small-aspect- 
ratio waves is derived. The model is shown to be completely inte- 
grable. The Lax pair, the first conserved quantities as well as the 
symmetries are exhibited. Theoretical and numerical studies reveal 
that it supports periodic progressive Stokes waves which peak and 
break in finite time. Comparison between the limiting wave solution 
of the asymptotic model and classical irrotational results is performed. 
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1 Introduction 



Wave propagation in an ideal incompressible fluid is a classical matter of in- 
vestigation in mathematical physics. Especially, surface gravity waves have 
been intensively studied and many model equations were introduced to han- 
dle this problem. Such studies are in a large part motivated by the fact 
that the initial three-dimensional water wave problem is not tractable an- 
alytically. Traditionally, the subject was studied in two almost separated 
domains: shallow water models and deep water models. 

The shallow water theory is widely known. The irrotational or rotational 
Euler equations have been approximated from several sides. An approach 
is based on the use of two parameters a and (3. a = a/h measures the 
amplitude a of the perturbation with respect to the depth h and /3 = h'^/X^ 
measures the depth with respect to the wavelength A. Assuming a and 
f3 smaller than one, a perturbative procedure is carried out. Then model 
equations are obtained by retaining only the lowest order terms in a and 
f3 pp. Another well-established approach for the derivation of shallow water 
model equations is the reductive perturbation method |2l |3l HI |5] based on the 
Gardner-Morikawa transformation [6]. It enables us to introduce slow space 
and time variables able to describe the effect of nonlinearity and dispersion 
asymptotically in space and time. An alternative method was introduced by 
Serre [7] and several years later by Su and Gardner [B] and Green and Naghdi 
[H [9] . It is based on an Ansatz which uses the shallow water limit of the exact 
linear solution of the rotational Euler equations. This way were obtained the 
Nonlinear Shallow Water [1] , the Boussinesq [TOl [11] , the Korteweg-de Vries 
(KdV), the modified KdV (mKdV) [121 US], the Kadomtsev-Petviashvili [H], 
the Benjamin-Bona-Mahony- Peregrine [IS], the Serre, the Green-Naghdi and 
more recently the Camassa-Holm fi7\ equations. All these model equations 
govern the asymptotic dynamics of wave-profiles of long-waves in shallow 
water. 

The deep water limit consists in considering an infinite depth (h = — oo) 
in the Euler equations. Then neither a nor /3 can be defined. The Gardner- 
Morikawa transformation, based on the dispersionless shallow water linear 
limit of the Euler equations, cannot be defined any more because the deep wa- 
ter limit is fully dispersive. Consequently, surface wave propagation in deep 
water is mainly concerned with the nonlinear modulation of wave trains. The 
most representative model equation is the ubiquitous NonLinear Schrodinger 
equation (NLS) [121 HSl US]- For a full account on modulation of short wave 
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trains in water of intermediate or great depth see Mei [20]. However, some 
model equations have been given for wave profile evolution in deep water. 
E.g., in reference [2T] is derived a finite-depth Boussinesq-type equation for 
the profile of irrotational free surface waves. The derivation is carried out 
via the theory of analytic functions followed by a perturbation procedure. 
Its deep water limit is exhibited as a function of integral operators. The 
Benjamin-Ono equation, first introduced by T. Benjamin [22] and H. Ono 
[23] , is worth being mentioned. It is a nonlinear partial integro-differential 
equation which describes one-dimensional internal wave profile evolution in 
deep water. In reference [23] an approach analogous the one we develop in 
the present paper was introduced. Other model equations whose solutions 
behave as deep water waves can be found in references [25] and [2S]- Their 
dispersion relations coincide exactly with that of water waves in infinitely 
deep water and the nonlinear terms are chosen in a ad hoc way in order to 
reproduce the Stokes limiting wave. 

The present work is the extension of a long sequel of studies carried out 
by the authors in various physical contexts, dealing with model equations 
for nonlinear and dispersive short wave dynamics ([2Z]-[3I])- Especially, the 
purpose of the paper is to study the dynamics of an elementary wave profile 
in deep water. Instead of looking for modulation dynamics of a wave train, 
we seek for the combined effects of dispersion and nonlinearity on a given 
Fourier component with wave vector k. Our approach is a generalization 
of the method firstly used by Serre for surface waves in shallow water. We 
assume an Ansatz based on the exact linear solution of the Euler system for 
deep water. Of course, as in the shallow water case, the Ansatz does not 
yield an exact solution of the full Euler system. It produces a depth- average 
model equation, which is expected to be valid close to the linear limit of the 
Euler system. 

The paper is organized as follows. In Section (|2]) we introduce the Eu- 
ler equations in deep water. In Section (|3]), after a brief summary of the 
"modus operandi" of the linear pattern Ansatz in the shallow water limit, 
we generalize it to the deep water case. In Section (jl]) we give the nondimen- 
sionalisation of the Euler equations and we introduce the deep water surface 
wave problem. In Section (|5]) we introduce the Serre approximation in the 
deep water context and we derive a nonlinear and dispersive system of equa- 
tions for the free surface elevation and for the free surface velocity. In Section 

we derive a small-aspect-ratio wave equation in deep water. In Section 
d?]) are exhibited the mathematical properties of the model: the Lax pair, the 
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relation with the integrable Bullough-Dodd model and the symmetries. In 
Section ([HD is studied the progressive periodic wave solution and the hmiting 
wave. The work is completed with a numerical study of the limiting wave. 
In the Appendix we give some detail on the derivation of the Bullough-Dodd 
equation and finally Section ([9]) draws the conclusions. 



2 The Euler equations in deep water 

We considerer the dimensional Euler equations for deep water. Let the par- 
ticles of the fluid be located relative to a fix rectangular Cartesian frame 
with origin O and axes {x,y,z), where Oz is the upward vertical direc- 
tion. We assume translational symmetry along y and we will only con- 
sider a sheet of fluid parallel to the xz plane. The velocity of the fluid is 
V{x, z, t) = {U{x, z, t), W{x, z, t)). The fluid sheet is moving on a bottom at 
z = —CO, and its upper free surface is located at z = ri{x,t). The continuity 
equation and the Newton equations (in the flow domain) read as 

U, + W, = 0, (1) 
pU + P: = 0, (2) 
pW + P: + gp = 0, (3) 

in which subscripts denote partial derivatives, and the dot ( ' ) denotes the 
material derivative defined by 

• dF OF ^BF ^^dF 

F = — = h U h W — . (4) 

dt dt dx dz 

The boundary conditions aX z = — oo and at 2; = T]{x,t) are 

W = for z ^ -00, (5) 
P* - Pq = for z = ri, (6) 
rjt-W + Urjx = for z = 7] , (7) 

where p is the constant density of the water, g the gravitation constant and 
Pq the atmospheric pressure. The solution of the problem consists in finding 
ri{x,t), V{x,z,t) and P{x,z,t). 
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3 A deep water linear pattern Ansatz 



Approximate solutions to the nonlinear surface water waves problem can be 
computed only if the z dependence of the velocity field V{x^ t) is known. 
The z dependence provides a coupling between the wave motion at the surface 
and the wave motion at various depths. In the linear case the z dependence 
is known, hence the solution is exactly known. Expansions in power series of 
of z as well as the theory of harmonic functions have been used to solve this 
problem in the nonlinear and potential Euler equations. Another approach 
to the problem is the use of an Ansatz based on the linear z dependence of 
V{x,z,t). It is the linear pattern Unsafe introduced in references [3 El El IS] , 
in the shallow water context, and equivalent to the widely known columnar 
hypothesis. It is thus worthy to begin with a brief description of how the 
Ansatz works in the shallow water case, before we extend it to the deep 
water context. 

The shallow water Ansatz. For surface progressive waves of wave vector k 
in water of depth h in the shallow limit kz ~ 0, the linear solution Ui{x, z, t) 
to the horizontal component U{x,z,t) of the velocity and the dispersion 
relation cu are 



Ui{x,z,t) 
uj{k) 



A ex^ i{kx — ut) 1 + Q{kz 



I- {khf + ... 
o 



kJgh 



(8) 
(9) 



where A is a constant, in complex representation so that the physical velocity 
is given by real part of Eq. ([S]). f/^ is a wave crest moving at the phase velocity 
c{k) given by 

c{k) = ^\l-l:{khy + ..]. (10) 

The linear pattern Ansatz assumes that the z dependence of U{x, z, t) is the 
same as the z dependence of Ui{x, z, t), and that its spatiotemporal behavior 
is given by an undetermined function u{x,t). So, Serre has assumed that 



U{x, z, t) = u{x, t), 



which corresponds to the extreme long wave limit in (|8]), i.e, U{x,z,t) does 
not depend on z. Since the phase velocity is invariant under k — )■ —k, i.e.. 



c{k) = c{—k), 



(12) 
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the Ansatz does not need to take into account the direction of propagation 
in u{x, t). 

The deep water Ansatz. We propose in the present paper an analogous 
procedure in the deep water case. Nevertheless two important differences 
arise. The first one is that in the deep water case, according to Eqs. ([IHT]), 
we have 

Ui{x,z,t) = Aexpi{kx — Qt) exp kz, (13) 

n{k) = v^, (14) 

c(fc) = (15) 

So, the phase velocity c{k) is no longer invariant under k — )■ —k and expres- 
sions (fT3l) . (IT^ and ffTSj) are strictly valid for wave crests moving to the right 
only . Consequently, the Ansatz would take into account the choice of the 
propagation direction made in deriving the linear solution. The second differ- 
ence is a very subtle one. Here we are interested in the effect of nonlinear ity 
and dispersion on a purely sinusoidal wave of wave vector k and frequency fl 
for large x and t. However, wave vectors and frequencies are no longer con- 
stant in dispersive system, and are generalized to the concepts of local wave 
vector k{x,t) and local frequency Q{x,t). There will then be a non-uniform 
local wave train (nearly sinusoidal) with k = k{x,t) and Q = fl{x,t). It can 
be shown that k{x, t) and Q{x, t) remain constant for an observer mov- 

ing at the group velocity Cg{k). Thus, to study the asymptotic behavior for 
large x and t of a wave with given wave vector k in deep water, the analysis 
must be carried out in the frame R(cg) traveling at the velocity 

c,{k) = lc{k) = (16) 



with respect to the frame R where the Euler equations ([T]17]) are written. 
Therefore, the linear pattern Ansatz we propose is 



U {x, z,t) = u (^k (^x — -ct^ , exp kz, c = c{k) = -y^ . (17) 

4 The deep water surface wave problem 

The first step is to nondimensionalise equations ([UlT]). Hence the original 
variables of space, time, velocity components and pressure will be normalized 
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with k, \/kg, yk/g and k/pg respectively, k is the wave vector of the wave 
to be described. Hence the Euler equations for — oo < z < 77 become 

U. + W, = 0, (18) 

u = -p: , (19) 

W = -PI - 1 , (20) 

and the boundary conditions are 

W = for 2 ^ -00, (21) 
P* = Po for z = r], (22) 
rjt = W — Urjx for z = rj. (23) 

We integrate the Euler equation in the depth of the fluid. From equation 
(fTSjl and boundary condition (1211) we get 

W{x,z,t) = - r U,{x,(,t)dC. (24) 

J —00 

This expression for W allows us to write the material derivatives U and W 

as 



U = Ut + UU, - / U^dC , (25) 

J —00 

w = - r u^tdc -u r u^^dc + r a^dc . (26) 

J —00 J —00 J— 00 

Using (1261) into (120|) . integrating with respect to z and making use of (122|1 . 
we obtain the pressure P*{x,z,t) as 

P* = r^-z-Tdclt U^tdC' + uf^ U^.dC-uJ^ UjA- 

Jz \j —00 J — 00 J —00 J 

+ Po (27) 



Note that expression fl27j) for the pressure diverges as z — ?■ —00 in Archime- 
dian way as it is expected to. The Newton equation f|T9l) reads, using fl27|) . 

as 

Ut + t/t/, - r U.JC' = - ^. + 1 r dC U^tdC" 

J~oo \Jz J —00 

+ f'udC u^^dc" 

J z J —00 

- f'u^dC f U^dC] . (28) 

J z J —00 ) X 



Now using Leibnitz's rule in the integrals in the right hand side we obtain 



Ut + UU,-uJ' U,dC = -v.+ TdCf U,,tdC" + V.r U^tdC 

J —OO Jz J~OD J —OO 

Jz J —OO J —OO 

- rU,.dC U^dC" 

J Z J — OO 

- v.uM r ujc , (29) 

^— OO 

with Uirj) = U{x, z = ri,t). Finally equation fl23|) expresses as 

Vt + U{r])v. + /' U^dC = . (30) 

J —OO 

If the z dependence of U{x, z, t) was known, the surface wave problem in 
deep water could be solved: Eq. ( 1271) gives the pressure P{x,z,t), while 
equations (l29l) (evaluated in z = r]) and (l30l) yield a system of coupled 
nonlinear equations for U{x,ri,t) and r]{x,t), which completes the solution 
of the problem. However this is not the case because U{x, z, t) is unknown 
versus z. 



5 A Serre-type system of equation 

The z-dependence of U {x, z, t) is the central issue making unsolvable the 
nonlinear surface water wave problem in deep water as in shallow one. Con- 
sequently we use the deep water Ansatz (ITTll . which reads in dimensionless 
variables as 

U{^,z,t) = u{^,t)expz, with ^ = x-h. (31) 

It implies the change in the differential operators: 

d 1 d d 
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Then we obtain from f l29p the expression 



Ut exp z — -u^ exp z + rj^ = — -u^^^exp^; + (-Ugfexpr])^ + -'u(M^^exp2r7)g 

— -uu^^^ exp 2z — -u^{u£^ exp 2r])^ 

1 1 , 

+ 2"€"«« ~ 2 ''"^^ ^''^ 

+ ^w^^exp^;. (33) 

This equation is now taken into account through an averaging over the full 
depth of the fluid, defined by the expression 

lim r{*}dC . (34) 

Z0-+-OO 1] — Zq Jzq 

Applying ( l3l|) to ( 133|1 we obtain 

{w^t exp(r7)}^ - i {u^^exp{ri)}^ + ^ {exp(2r7) [mm^^ - m|] }^ - = 0. (35) 

Finally, substituting (EI]) into fl5U]) we get 

Vt - + {u exp r])^ = . (36) 

Note that we did not assume that either ri{C,,t) or u{C,,t) vanish as — )■ oo. 

The system fl35|) and f p6|) for r7(,^, t) and u(^, t) is a z-average of the Euler 
equations resulting from the flow hypothesis fl^ and the subsequent average 
in the depth f l5^ . It is a deep water analogous of the Serre (or Green-Naghdi 
or Su-Gardner) system. It differs from the result obtained in [2l] . The latter 
is a weighted z-integration of the Euler equations. The crucial step there was 
the use of a precise weight in order to regularize the Archimedian divergence 
term present in the pressure. In turn this produces a weight-dependent model 
the degeneration of which is eliminated by the requirement that the linear 
limit possesses phase and group velocities equal to those of the deep water 
system. The present derivation does not need either the weighted average 
procedure or the ad doc hypothesis of equality between the linear limits. 
Both of them, z-average and weighted z-integration, coincide in the classical 
shallow water case. 



9 



Assuming that ?7(^, t), t) vanish for ^ — )• oo, and setting Oq = dt — \d^^ 
system ( |35ll and (1361) can be reduced to the integro-differential model for 

1 2 3 /•« 1 

Vee + V- 7;Ve - V( / Veed^ - T^V^e / Ved^ + ^ 

^ J —OO L J —oo L 

(37) 

This equation for r]{^,t) is close to classical formulations of model equations 
for surface waves in deep water, which are integro-differential models. 




6 An small-aspect-ratio surface wave model 

Up to this point the asymptotic behavior in ^ and t of the surface wave 
profile was not considered. This can be performed by taking advantage of 
some small parameter present in the system. Therefore, let us rewrite the 
nondimensional function t] as 

r]{^,t) = kaH{^,t), (38) 

where the nondimensional function if(^, t) represents the surface deformation 
with amplitude a measured from the initial mean water level z = 0, and the 
product 

ka = e , (39) 

is the steepened coefficient or aspect ratio, which measures how much the 
wave profile is peaked (e > 1) or flat (e < 1). The equations fl35l) and fl36l) 
are written explicitly in terms of e and H as 

exp(eif)| + |iexp(2eif) [mm^^ - eH^ = 0,(40) 

eHt-^eH^ + [uexp{eH)]^ = 0.(41) 

In fact, pUj) and fHTl) are valid for any orders in e, but they are dispersive 
and strongly nonlinear, so untractable analytically. Nevertheless, a pertur- 
bative theory can be carried out in terms of e by considering flat propagating 
waves, i.e., e < 1. The linearized system yields the nondimensional form of 
the dispersion relation 

n = ^. (42) 
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Nonlinearity causes deviations from the latter (Stokes' hypothesis) which can 
be taken into account through 



= - + e + 



(43) 



Hence the phase becomes 



^ - nt = ^ - -t - et - eH + ... , 

and we can define new variables y, r, z/, as 

. 1 2 
y = ^--t, T = et, iy = et,.... 



(44) 



(45) 



The function if(^,t) must then be regarded as a function of the new inde- 
pendent variables y, r, z/, ... jlH |44], i.e., as H{y,T,h', ...,). On the other 
hand, the derivation operators with respect to x and t become 



d _ d 

dy' 



d_ 

di 



1 d 



d 

+ e— + 



2 dy dr 



d_ 
du 



(46) 



Since exp(eif) = 1 + eH + O(e^), we can retain in equations ( l40l) and (l4T]) . 
after the leading order e°, both orders e° and e, and so on. Thus, Eq. ( l4Tll 
reads as 

^9j, + e9, + 0(e2)]eif-i9j,eif+ [M(l + eif + 0(e2))] =0. (47) 



This equation gives u in terms of derivatives and antiderivatives of if, as 



u = eH 



Hrdy' 



O e 



(48) 



where was not assumed that u ot H (and derivatives) vanish as y — )■ — oo. 
Instead of such assumption, the above equation implies the relation 



u{—oo, t) = eH{—oo, r) 1 — eH{—oo, t) + O (e 



(49) 



Reporting the expression of u into (jlQj) and retaining terms in e and e^, we 
obtain the equation 



- 2e^Hryy + eHy + eifj^^^ - -e^{HHy)yy = 



(50) 
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which in terms of the dimensionless variables and ri{^,t) reads as 




3 

'^Vi^t = V^- ^ivVi)ii > (51) 
or in terms of variables with dimension: 

3 

= k'^V( - i^Kmd^i ■ (52) 

This equation is the main result of our work. It describes the asymptotic 
nonlinear and dispersive evolution of small-aspect-ratio waves of a Fourier 
wave vector k in deep water. It has /^-dependent coefficients. Equation (1521) 
can be considered as belonging to both of the two categories: that of KdV 
models (KdV, modified KdV, Benjamin-Bona-Mahony-Peregrine, Camassa- 
Holm, etc.) describing evolutions of wave profiles and that of NLS-type 
equations (modified NLS [13], Davey-Stewartson [35], etc.) describing mod- 
ulation of wave profiles and having fc-dependent coefficients. 



7 Mathematical properties 
7.1 Lax pair 

In this section we assume that rj and its derivatives vanish as — ?■ icxo. 
Equation (ISTl) has the dispersion relation 

= 1 , (53) 

as required, because ( 15T|) was derived in a frame moving at the group velocity 
associated with the Euler equations in deep water. Making use of the scalings 

^\,(i.t) , 2|^|, (54) 

mt = Vi- (w?)ss, (55) 



Eq. (l5T!) becomes 
and integrating once yields 



V^t = V- ivVi)^ ■ (56) 
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Let us consider the function F defined by 

F' = l- 3r]^^ , (57) 

or in terms of (p = rj^^, by 

F^ = 1- 30. (58) 

One of the most remarkable properties of F is that it allows us to build 
non-trivial conserved quantities. The two first ones are 

Ft =- ivF)^, 
(20|F-^)^ = [(l + 0)F-i-2#|F-^]^. (59) 

The Lax pair also can be constructed by means of F and reads 

M ^ (60) 

where di, (72 and cfs are the Pauli matrices, A the spectral parameter and Eq. 
f lSB]) is obtained as the classical Lax equation Lt = [L, M]. 

7.2 Reduction to the Bullough-Dodd equation 

Through the change of variables from to p and the change of functions from 
ri{^,t) to r{p,t) defined by 

P = f^Fd^, (61) 



1 - 3?^^^ = exp(r) , (62) 
it is found that r{p,t) satisfies the Bullough-Dodd equation (see Appendix) 

rpt = ^ jexp Qr^ - exp | • (63) 

The Bullough-Dodd equation was introduced in [35], in which some non- 
trivial conserved densities were also exhibited. Its Lax pair and complete 
integrability were shown in ^7]. The change of variable and of function 
fl6T]) and fl62|) leading from fl56|) to fl63l) were introduced in the context of 
a systematic study of short-wave dynamics in long-wave model equations in 
references |1S1 [HI EDI EI] ■ 
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7.3 Symmetries 

Equation fl55|) is invariant under the discrete transformation 

v{U)-^v{-^,-t) , (64) 
and can be written as the conservation law 

dtiv^() = d^{7] - {r]r]^)^) . (65) 
It is invariant under the Gahlean group of transformations 

t' = t, x' = ^ + Vt, r]{^,t) = 7]o + S{^,t), (66) 

in which V and rjQ are constants with rjo = —V. In this case indeed, Eq. f l55|) 
transforms into 

Sx'x't' = Sx' — {SSx')x'x' ■ (67) 

Consequently, (1^ must be seen as a member of a family of equations 
parametrized by the speed V of the Galilean group of transformations fl66p . 
This invariance comes from the invariance under fl66l) of the Euler equations. 
Last but not the least, the equation (|55|1 - or (|56|1 - has the Lorentz invariance 

^ ^ <, t ^ -, f^'^V, (68) 

with K, a real arbitrary parameter, which is undoubtedly related to the inte- 
grability of the associated Bullough-Dodd equation. 

8 Periodic wave and limiting wave 

8.1 Harmonic solution and perturbative approach 

A periodic wave solution of f lS^ is predicted by linear analysis. It is sinu- 
soidal, of small amplitude and travels with phase velocity {l/2)^g/k. The 
question is: up to what values of its amplitude and velocity does the peri- 
odic wave exist? This is still possible for wave amplitudes large enough so 
that higher terms in the perturbation series can no longer be neglected in 
Eq. ( 152|) . Then the wave is not exactly sinusoidal, its velocity is not exactly 
^^g/k, and the limiting wave has the greatest height before breaking. In 
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order to find this steep rotational Stokes wave, let us start with equation fl52l) 
in the frame R with coordinates {x,t), which reads as 

[k 2 9 3 

~Vxxt — k Tjx — Tjxxx ~ '^kr]x1]xx ~ '^krjTjxxx ■ (^9) 

The radius of curvature is defined by 

Vxx 

Hence Eq. fl69|) . within the allowed order, can be written as 

jk Rt 2 , -^a; QjVx , ^, Rx /-.x 

We assume a progressive periodic wave in the variable z = x — ct, hence we 
have 

'2c^Rz = k^R-" + R,- ^kv.R + hvRz ■ (72) 
Let 7] be the height of the wave, i.e., 

V Vcrest Vtroughy 

(73) 



max 



and assume that the wave approaches the limiting wave of height rj = t], 
and that this value is reached for z = zq; i.e, ri{zo) = rjmax with zq = xq — c^to 
where cl is the phase velocity at the limit. At the limit, on the one hand, the 
radius of curvature -R(^o) is zero [Ml [51], and on the other hand its derivative 
Rz{zq) is not well defined. Nevertheless, equation (172]) can be satisfied if 

[k 3 

'2CL]l- = l + -kr]max- (74) 

To determine the limit krjmax and cl, one further equation is required. This 
equation will come from the Stokes series for the periodic solution of small 
amplitude of Eq. ( 169]) . Hence we consider a periodic solution 

r]{z) = 6{7]o{z) + 6r]i{z) + 6%{z) + ...} , (75) 

with S a small amplitude defined as 

^ = ^{Vcrest - Vtrough) = (76) 
15 



To avoid secular terms, the velocity c must be also expanded as 

c= cq + 5ci + 5'^C2 + .... (77) 

We write Eq. (1^ in terms of report expansions (1751) and (177|) into it, and 
then isolate order by order the coefficients of each power of 5. We obtain this 
way a series of equations, whose leading order 5 is 

(^1 - 2co^^ r?o,3. - k'^r^o,. = ■ (78) 
An even solution of fl78l) is 



?7o(-2) = cos(A;z) with '^o ~ y C^^) 
Using the solution (179|) . the equation at order 5^ is 

^?i,32 + = -2ciyf P sin(fc2) + 3fc^ sin(2fcz) . (80) 

The trouble is that the term in sin(A;2;) resonates with the operator on the 
left-hand-side. This is a secular term. It can be eliminated by choosing 
ci = 0. Thus the solution at this order is 

k 

T]^(z) = -cos{2kz) with ci = . (81) 

At order 6^ we find 

2c2W -A;=^ - — j sin(A;2;) - sm{3kz) , (82) 
and the solution free of secularity is 

miz) = '^^^cos{3kz) with C2 = ^y^. (83) 

Thus, the solution t]{z) at order 6^ is 

k 27fc^ 

ri{z) = 6 cos{kz) + 6'^— cos{2kz) + 6^ cos{3kz) + 
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Now, when the solution approaches the hmiting wave solution we must 
have 

6 — y ^"^^ and c — > cl, (85) 

hence 

This is the second equation we were looking for. From ( 1741) and ( l86l) we 
obtain 

kVmax = 0.697, = 1.022 . (87) 

In the irrotational case the classical result found a long time ago by J. H. 
Michell [55] for krjmax and recent computations for cl [56] yield the values 

kvmax = 0.892, cJ- = 1.092. (88) 

V 9 

Comparing (IHTl) and (!88l) we find the percentage errors (100% times the 
relative errors) 

{k7]^ax) error = -21%, (^^V-j ="6,4% (89) 

\ ) 9 / ^rror 

Deviations of krjmax and cl, calculated from (!69ll in relation of ( l88l) are due 
to three mains factors: 

• the results (!87|) come from a computation at second order in 6 only, 
while (18 8 p were obtained using theories at least at fifth order in S, 

• the classical results were obtained from the complete Euler equations, 

• our fluid is rotational. 



8.2 Exact analytical approach to the limiting wave 

We start with the equation (15 ip in the referential R relative to the variables 
(^,t), i.e., 

2r]^t = V-m-liv\i- (90) 
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We look for a traveling wave with speed c, 



Tj = rj{z) , z = k{C, — ct) 

so that fl9UI) can be written as 

■3 



-7] +ari 



(91) 



(92) 



where we have set a = 1 — 2c. Now multiplying both sides by 

-3 



-7] +ar] 



and integrating once, we obtain 



-T] +ar] 



if + arf + K, 



K being some integration constant. And hence 

?>r] + 2a 



dz 



■drj. 



"2a/?7^ + aT]'^ + K 
Setting 7] = aY, z = ZJ\a\ and K = \a\^p transforms (!94l) into 



(93) 



(94) 



dZ 



3y + 2 



2Jp + £ (r3 + y2 



--dY. 



(95) 



e being the sign of a. Then we look for some regular periodic solution. It 
attains its minimal value Yi for some Z, say Z = Zi, and then will grow 
up to some maximal value reached aX Z = Z2 for the first time after Zi. 
dY/dZ must be well-defined, real and positive for all Z between Zi and Z2 
and all Y between Yi and Y2. Further dY/dZ must be zero for Z = Zi and 
Z2. The conditions Zi < Z^ and Yx < Y2 must also be satisfied. 

fiY) = y3 + F2 presents two local extrema, /(-2/3) = 4/27 and /(O) = 
0. Hence p + sf{Y) has more than one zero if p is between and — 64:121. 
Since dZ/dY is zero for Y = —2/3, and dY/dZ must remain finite one the 
whole interval, the value F = — 2/3 is excluded, and Yi and Y2 must be the 
two largest zeros of p + ef(Y), with Yi < < 1^2- For e = +1, p + sf(Y) is 
negative in this interval. Hence, since its square root must be defined, the 
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Figure 1: The normalized amplitude Y2 — Yi of the Stokes wave vs the 
parameter p. 



case £ = 1 is excluded. The period of this solution, which is the wavelength 
of the Stokes wave in normalized units, is 

The amplitude Y2 — Yi is given by the solution of the cubic polynomial 
equation + — p = 0, it is plotted vs p in Fig. [T] The integral ( PB]) 
can be computed numerically, which yields the normalized wavelength vs the 
integration parameter p as shown on Fig. |2j Finally the wavelength is shown 
versus the amplitude in Fig. |3l The maximum amplitude it attained with the 
minimal wavelength for p = 4/27. Then the quantity p — Y"^ — Y^ factorizes 
explicitly, and after simplification we get 



dZ 3\/3 1 



(97) 



dY 2 V1-3F 
which is straightforwardly integrated, and then inverted, to yield 

y = - - — , (98) 

3 9' ^ ' 

in which we have set arbitrarily Zi = 3 (translation invariance). Formula 
fl98|l is completed using parity and periodicity to 
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Figure 2: The normalized wavelength — of the Stokes wave vs the 
parameter p. 




Figure 3: The normahzed wavelength Z2 — Zx of the Stokes wave vs its 
normalized amplitude Y2 — Yi. 
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E denoting the integer part. Coming back to the original variables, we obtain 
the expression of the limiting Stokes wave, as 



^ 36 



2z (z 1 



(100) 



in which A = 6a/— a is the wavelength (it has been seen that e = sign (a 
— 1). The velocity is 

c - ^ " - V 



2 2 2 72 

The total amplitude is —a = A^/36. 

Coming back to the dimensioned variables, we obtain 



(101) 



kX 



l2 



V 



36 



(102) 



in which X' = X/k is the dimensioned wavelength, and 



c 



X 



72 



-Vat. 



(103) 



We need to identify k. The exact linear wave must be periodic with wave- 
length Xi = 2TT/k. Assuming the same value of k for the same wavelength in 
the nonlinear case (A' = A;), we get (dropping the prime) 



27rA 



(104) 



and ( = X — cit, with 



(105) 



The limiting Stokes waves, in the dimensionless form fl99|) . is plotted on 
Fig. m Its dimensioned form ( 1104p is plotted on Fig. [5] for A = 1. Recall 
that Y/i] is negative! Now the percentage errors using fll04p and fllOSp are 



{krj, 



max J error 



+22%, 



-4%. 



(106) 
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Figure 5: The limiting Stokes waves in dimensioned form, for A = 1. 
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8.3 The evolution of the Umiting Stokes wave 



The evolution of the hmiting Stokes wave is checked numerically. We use the 
Eq. fl90l) . written in frame traveling at velocity 1, as 

2v^t = V-l{v')i^- (107) 

Eq. (I107p is integrated with respect to x, then we apply the Fourier transform 
defined as 

y{^,t) = Jy{a,t)e'''"'^d^, (108) 

to get 

1 3 - 

Vt = -r- — ^ ~ -riTxcrr]'^ ■ (109) 
4«7r(T 4 

Then the time evolution of t) is computed by means of a standard 4th or- 
der Runge-Kutta algorithm. The Fourier transforms are computed using a 
standard fast Fourier transform algorithm. It must be noticed that the term 
representing the antiderivative in Eq. fllOQp is not defined for a = 0. The 
corresponding term in the discrete scheme is set to zero, which assumes a 
zero mean value. The mean value of the limiting Stokes wave (11041) is easily 
computed, it is zero as required. Further, the fact that (11041) is not continu- 
ously derivable induces a high-frequency numerical instability. The latter is 
removed using a spectral filter in the numerical scheme, which attenuates the 
highest frequencies. It is checked that this filter does not affect the spectrum 
of the limiting Stokes wave itself. The result of the computation is plotted 
on Fig. [61 



9 Conclusions and final comments 

We have obtained a deep water equivalent of the widely known Serre shallow 
water system of equations. This was done via an appropriated Ansatz for the 
vertical coordinate dependence of the velocities. This procedure allowed us 
for a description of the nonlinear and dispersive dynamics of the wave profile 
itself, avoiding the classical route of wave modulation of a short wave train. 
From the system, we were able to derive an asymptotic model equation (in 
times and space) for short waves. The new model is the short wave anal- 
ogous of KdV in surface waves and, as KdV, shown complete integrability. 
The progressive Stokes wave was studied analytically and numerically. The 
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Figure 6: The evolution of the hmiting Stokes wave, as computed numerically. 

results for the limiting wave height and the limiting phase velocity, within 
the degree of approximation of the model, are in accordance with classical 
results. Important future works and open issues to be study are : a) the anal- 
ysis of the Stokes limiting wave of the system for the free surface elevation 
and the free surface velocity, b) inclusion of surface tension, small viscosity 
as well as higher spatial dimension, c) the achievement of asymptotic sys- 
tems of higher degree than two in the aspect-ratio parameter, d) is the new 
asymptotic model the first term of a hierarchy of integrable models, as KdV 
is? e) last, but not the least, an important open issue is: can we go fur- 
ther in the unexpected bridge established by the system between nonlinear 
and dispersive short wave dynamics in fluid mechanical and two-dimensional 
integrable relativistic field theory? 
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Appendix A. Reduction to the Bullough-Dodd 
equation 

The change of variables from ^ to p and the change of functions from ?7(^, t) 
to r{p, t) 

P = t Fdi, (110) 



F = (1-3^5^)3, (111) 
l-3r/55(e,t) = exp(r(y,t)), (112) 

together with the conservation of F 

Ft = -ivFh (113) 
gives for any function M(^, t) the relations 

dp ' 

^ = + (115) 

at op at 

From (11151) follows that 

8^mA ^'^(^^^) ^ dM{p,t) , ^ d^M{p,t) 

^ ^ '"^'^^^ + ^^^^^ 

We take the derivative of the evolution equation fl56l) with respect to ^, once 
and then twice. Using the definition of F, it yields 

+ = Vl;F^ (117) 

+ = mF^ - ^VtiFri^ip. (118) 

Equations (11141) and (11151) for M = r]^^ are reported into (I117p . which gives 

1 / N 

^C = ;^?/«t- (119) 

Then using (I116p for M = rj^^, and (IllSp and (I119p . brings to the important 
relation 

V^^ptF^ = F\^ - 3r]^^tv^^p. (120) 
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The change of function fl61f62p and the expressions of F, Fp and Ft give 

fl211 



p 

3. 



1 J F\i:pt + 3r]^^tViiip 
^* ~ 3 F6 



Finally using f ll20p and the expression of F we obtain the Bullough-Dodd 
equation 
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